Universal statistics of the critical depinning force of elastic systems in random media 
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We study the rescaled probability distribution of the critical depinning force of an elastic system 
in a random medium. We put in evidence the underlying connection between the critical properties 
of the depinning transition and the extreme value statistics of correlated variables. The distribution 
is Gaussian for all periodic systems, while in the case of random manifolds there exists a family of 
universal functions ranging from the Gaussian to the Gumbel distribution. Both of these scenarios 
are a priori experimentally accessible in finite, macroscopic, disordered elastic systems. 
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The study of disordered elastic objects sheds light on 
the physics of numerous experimental systems, including 
domain walls in magnetic [jj] or ferroelectric j^j] materi- 
als, contact lines of liquid menisci on a rough substrate 
and propagating cracks in solids ^ . All these are ac- 
curately modeled as elastic manifolds in the presence of 
randomness. A second class of disordered elastic systems 
is given by periodic structures such as charge density 
waves i3] , vortex flux lines in type-II superconductors 6] 
and Wigner crystals It has been proven that the 
periodic nature of these structures carries to the same 
universal behavior as elastic interfaces in the presence of 
periodic disorder [1,0. 

A continuing challenge is to understand the response 
of these systems to an applied external drive per unit 
surface, /. Two regimes are present at zero temperature: 
whenever / is smaller than certain critical threshold fc 
the interface is pinned, but beyond this poi nt ( f > the 
system undergoes a 

and continues to move with a non-zero average veloc- 
ity. Throughout all of the low temperature regime the 
dynamics is strongly influenced by the presence of this 
threshold. It is thus interesting to characterize fc in 
detail, and in particular, to study its sample-dependent 
fluctuations. We remark that the critical force is actually 
the biggest one that still pins the system. New insights 
shall then be expected from exploring the eventual con- 
nections with the theory of extreme value statistics. 

In a way reminiscent to the central limit theorem for 
the sum of uncorrelated variables, there exists a funda- 
mental result for the statistics of extremes [3 [13. It 
states that only three limits for the maxima of inde- 
pendent identically distributed (iid) variables are pos- 
sible: the WeibuU family (for upper bounded variables), 
the Frechet family (for power-law-tailed variables) or the 
Gumbel distribution (for exponential-tailed variables). 
However, due to the complicated dynamics of elastic sys- 
tems, we expect that strong correlations should play a 
crucial role in determining the limiting distribution of its 
thresholds. Indeed, the depinning transition is an exam- 
ple of a critical phenomenon with well-determined uni- 
versal exponents, which implies strong, long ranged, well- 



characterized correlations. The few results known related 
to extreme values of correlated variables concern short 
range or weak correlations J^, JJ,, J^J . The opposite 
case has only been started to be discussed 0, [13, [^ ■ 
Here, we identify and study an example of extremes of 
strongly correlated variables with clear physical interest. 

To pose the problem, let us consider a fixed disorder 
realization on a finite system. To every configuration of 
the elastic interface, {a}, we associate a depinning force, 
fjl°'\ given by the smallest non-negative force so that at 
least one point of the interface has instant non-negative 
velocity . The critical force of the realization fol- 

lows naturally: 



= max {/]">} 

{a} 



(1) 



In this letter we will concentrate on the statistics of 
We shall identify the universal aspects of its distribution 
and stress the important differences between the periodic 
and random manifold cases. We shall find that whilst for 
periodic systems the distribution is always a Gaussian, 
for random manifolds it belongs to a family that inter- 
polates between the Gaussian and the Gumbel "s^]. 

The focus of our attention will be a d-dimensional in- 
terface propagating transversally on a ((i+l)-dimensional 
space; for clarity we describe now the case of a line 
(d = 1). At any given time, the line defines a single- 
valued function h(x), where x is the longitudinal coor- 
dinate. Its zero temperature overdamped dynamics is 
governed by the following equation of motion. 



dth{x,t) 



dE{{h,x}) 
dh(x, t) 



cV^h+f + riix,hix,t)) . (2) 



The functional E{{h,x}) represents the total energy, in- 
cluding the harmonic elastic energy and the potential en- 
ergy terms due to the drive / and the short-range corre- 
lated disorder force rj^x, h). This is a non- linear equation 
whose analytic solution is unknown. As alluded to in 
the beginning, properties of the disorder distinguish two 
different cases. If ri{x, h) is a periodic function, Eq. (O 
describes a periodic system, whereas the random mani- 
folds are characterized by a non-periodic disorder. 
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Our approach to studying Eq. (0) will be numerical. 
In particular, for each disorder realization, we find the 
critical force using an exact algorithm Q that lo- 
cates the critical line h^{x) {i.e. the ultimate pinned 
configuration). Fig. schematizes the model: we keep 
h as a continuous variable and discretize the coordinate 
X ^ i] periodic boundary conditions are used so that h 
has period M and i has period L |3 Ell ■ For a periodic 
disorder M is kept constant as L goes to infinity. In this 
case the interface winds around the system and feels the 
disorder periodicity. For a random manifold both L and 
M go to infinity. The most physical scaling corresponds 
to M ~ L^, where C is the roughness exponent at the de- 
pinning transition, and allows the transversal dimension 
to scale with the lateral extension of the interface 0, ^3 ■ 

A general theorem for disorder-controlled transitions, 
shows that on a finite- volume any divergent cor- 

relation length scales with an exponent 0, Ull Eg 

..s > ^ . (3) 

The identification between vps a-nd the critical correla- 
tion length exponent is a very delicate problem that 
has been discussed in the literature In any case, for 
elastic systems it was determined that ^ |ll|, [l^, ll3 1 

al^W^W -L-^ , (4) 

where the over-bar indicates disorder average. 

Let us start discussing the periodic case, where by def- 
inition = 0. We have performed multiple simulations 




M = const. M = kL'^ 



FIG. 1: Elastic lines in a disordered medium. The left schema 
illustrates a periodic system whereas the right one shows a 
random manifold. 
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FIG. 2: Rescaled distributions for a periodic system. The 
squares correspond to the results for a periodic system (M = 
8) with L — 64: (open) and L = 1024 (filled); undisplayed error 
bars are smaller than the symbol size. The lines are the Nor- 
mal (solid) and Gumbel (dashed) distributions. Convergence 
with L towards the Gaussian limit is evident. The values for 
the skewness (71 < 0.03) and kurtosis (72 < 0.01) are consis- 
tent with this limit. Inset: covariance of the 77(i,/ic)'s (filled) 
compared with that of a set of random numbers (open). 

for d = 1, determining the critical force for different dis- 
order realizations. We verify that the scaling relation 
given in Eq. is obeyed with vps ~ 2 (in agreement 
with the results of previous simulations |^). In order to 
highlight the universal behavior in the thermodynamic 
limit, it is convenient to rescale the forces and compare 
always zero-mean unit-variance distributions, p(/c) with 

= (/_r — fl)/af^. In Fig. |21we show that the limit dis- 
tribution of p(/c) is Gaussian. Finite size effects are very 
small and the limiting behavior is robust against changes 
in the nature of the disorder, the period M , the boundary 
conditions or any microscopic parameter like the elastic 
constant (in all figures, the number of samples ranges 
in 10^-10^). Simulations done for d = 2 also confirm 
this scenario. The generic nature of this result suggests 
that, in all periodic systems, the limiting distribution of 
thresholds (critical field, current, etc.) is Gaussian. 

To understand this result, we study the correlation of 
the critical disorder- forces, ri{i,h^), defined as the pin- 
ning forces on the different points of the critical line. 
Since the configuration is static, the critical force is ob- 
tained as minus the average pinning force. The degree 
of correlation among the ri{i,h^ys plays thus a defining 
role in determining the final thermodynamic distribution 
of critical forces. Evidently, if they were uncorrelated, 
the distribution would surely be Gaussian, on the gen- 
eral grounds of the central limit theorem. We studied 
therefore their covariance, defined as 



cov(^) = E + ^' + fc){v{h K) + fl) ■ (5) 
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FIG. 3: Rescaled distributions for a random manifold. The 
circles correspond to the results for a scaling with AI — L'' for 
L — 64 (open) and L — 768 (filled). The lines and error bars 
are as in Fig. |5| Convergence with L towards a non-trivial 
limit is evident. The extrapolated values for the skewness 
and kurtosis are 71 0.57 ± 0.01 and 72 — » 0.63 ± 0.03. 
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FIG. 4: Cumulative distributions of critical forces of lines 
with L = 256. Different solid lines correspond to different 
A/'s as indicated. The average width is to ~ 100 {kw ~ 0.1). 
The rescaling according to Eq. (|SJ of the curves for M = 
64, 128, 256 to M = 2165 (A: = 2), is shown with dashed lines. 



We found that, in the case of fixed finite M, the cor- 
relation decreases very fast (at least exponentially) and 
the variables quickly become uncorrelated (see the in- 
set in Fig. 121), explaining the observed convergence of 
p(/c) towards a Gaussian distribution. The ri{i, h^Ys self- 
average in the thermodynamic limit, in turn implying 
that fps = 2/(i, saturating the bound of Eq. (Q. 

We direct now our attention to the case of random 
manifolds. We find i'fs = 1-33 ± 0.01, verifying for d — 1 
the scaling relation vya = = (2 — C)~^ [lOillllll^i with 
C ~ 1.26 ± 0.01 mliil. We stress that this relation is 
valid, as will be clearly shown later, only in correspon- 
dence to the physical scaling M ^ (cf. Ref. H^). In 
this case the distribution is non-trivial; our results for 
a d = 1 random manifold are presented in Fig. 13 We 
find a distribution that converges as i ^ 00 to a uni- 
versal curve that clearly deviates from the Gaussian in 
the direction of the Gumbel. Previous studies on cellular 
automaton models have found that distributions of near- 
threshold depinning forces are also universal [30lf3ll|. 

The proper scaling for the random manifold case 
should be defined with care. We consider M — kL'^ with 
k a constant independent of L. For each fixed fc, we find 
that the distribution of critical forces has a non-trivial 
thermodynamic limit. In particular, for fc — > this limit 
tends to a Gaussian, as seen in the periodic case, and in 
the opposite limit {k — > cxd) the distribution approaches 
Gumbel's. We verified this behavior for d = 1 and results 
for d = 2 confirm the same scenario. The most natu- 
ral definition of a random manifold system corresponds 
to considering values around certain, model dependent, 
kyj = w/ L'' (where w is the average width of the elastic 
line 112 )■ The ratio k/k^ then parameterizes a universal 



family of functions that crosses over from the Gaussian 
to the Gumbel distribution. This set of extremal distri- 
butions is universal in the sense of being independent of 
all microscopic parameters of the model and fixed solely 
by d and the ratio k/k^- 

To understand the Gumbel limit, let us start from the 
problem of iid variables. Given a set of n variables with 
Zn its maximum, the associated cumulative probability 
distribution, Hn{z) — P{z > Zn), obeys: 

Hniz) = [iJ™(z)]" . (6) 

This is the crucial property that allows to prove the the- 
orem on the statistics of extreme values mentioned at 
the beginning. We can now go back to the case of ran- 
dom manifolds. When the aspect ratio of the samples 
becomes very wide, k ^ k^, the system can be thought 
to be split into k/kw approximately independent regions. 
This is because the configurations with non-zero depin- 
ning force proliferate and their correlation decreases. In 
this case Eq. ® is asymptotically verified and the results 
of the theorem are still valid as shown by the statistics 
of extremes of staUonary sequences U M "^er 
to illustrate this point, we study explicitly the cumula- 
tive probability of critical forces, HM{fc)- We present 
in Fig. 21 curves corresponding to different values of M 
for L large enough that the finite size effects are indis- 
tinguishable in the scale of the figure. We used Eq. © 
to scale the distributions towards one corresponding to a 
much larger M. Let us separate our discussion in three 
regimes: k < k^, k > k^j and k ^ kw The figure shows 
that whereas in the first case the scaling of Eq. © is 
violated, in the other two it is obeyed with increasingly 
better accuracy as k grows. This demonstrates that, in 
the third case, the critical force can be thought of as the 



4 



maximum among forces sampled from a distribution like 
that in Fig. O Since the tail behavior of the latter is 
bounded between Gumbel and Gaussian (both belonging 
to the domain of attraction of the Gumbel distribution 
(isf). it follows that in the k ^ oo limit the distribution 
should be Gumbel, as opposed to Frechet or WeibuU. 

Let us discuss briefly the implications for other scal- 
ings of M. At this stage it should be clear that, for 
any other scaling, the distribution converges to a non- 
parametric thermodynamic limit. Consider, for the sake 
of concreteness, M ^ L'^ : on the one hand, for (' > C, 
the convergence is towards the Gumbel distribution; on 
the other hand, for < the limit is the Gaussian. 

As we have seen, the study of the distributions of 
maximum thresholds unravels a very rich and complex 
scenario. For all periodic systems the distribution is a 
Gaussian, -a non-trivial result from the viewpoint of the 
theory of extreme value statistics. This result rests on 
the asymptotic validity of the central limit for the crit- 
ical depinning forces and implies i^fs = 2/d. For ran- 
dom manifolds the situation is more intricate and there 
exists a universal family that interpolates continuously 
from the Gaussian to the Gumbel distribution. To our 
knowledge, this is the first time a set of distributions of 
extremes with these characteristics was identified, (there 
is another recent example of a family that interpolates 
between these two limits, but with no interpretation in 
terms of extremes 'SS^ ) . The criticality of the system is 
crucial to obtain such varied behaviors. Indeed, it is the 
infinite range of the correlations that build up along the 
critical manifold that is responsible for the strong devia- 
tions from the standard theory for iid variables. 

It would be extremely interesting to verify these sce- 
narios experimentally. We believe, good candidate mea- 
surements in periodic systems would be those of super- 
conducting critical currents, or charge density waves. For 
propagating interfaces the usual set-ups involve M 3I> w. 
To build a distribution in this case one should be able to 
measure the critical thresholds for a given tiling of the 
sample, then compare distributions among different tiling 
widths. On a different line, we already started to look 
into the analytical characterization of the family of dis- 
tributions that we found for the random manifold case. 
This might provide deep insights into the theory of ex- 
tremes of correlated variables that is still in its infancy. 
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